Structure of boson systems beyond the mean-field 
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Abstract. We investigate systems of identical bosons with the focus on two-body 
| correlations. We use the hyperspherical adiabatic method and a decomposition of 

the wave function in two-body amplitudes. An analytic parametrization is used for 
the adiabatic effective radial potential. We discuss the structure of a condensate for 
arbitrary scattering length. Stability and time scales for various decay processes are 
i estimated. The previously predicted Efimov-like states are found to be very narrow. 

We discuss the validity conditions and formal connections between the zero- and finite- 
\^ ' range mean-field approximations, Faddeev-Yakubovskh formulation, Jastrow ansatz, 

I/"") . and the present method. We compare numerical results from present work with mean- 

field calculations and discuss qualitatively the connection with measurements 
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1. Introduction 



The interest in dilute Bose gases has been growing since the experimental realisation 
of the phenomenon of Bose-Einstein condensation (BEC) plEUHj. Excellent reviews of 
the world of BEC are given in recently published monographs [3J E] . The theoretical 
interest in BEC goes more than fifty years back and is widely based on the mean-field 
formulation. The usual measure of the validity of the mean-field is that n|a s | 3 <C 1, 
where n is the density and a s is the two-body s-wave scattering length [3]. Then the 
particles are not too close to each other and correlations are expected to be negligibly 
small. The importance of correlations must increase with the density of the system and 
the mean-field method sooner or later becomes inadequate. 

A Feshbach resonance is routinely used to create Bose-Einstein condensed systems, 
where the effective interaction corresponds to a large scattering length a s 6\. Then 
stronger correlated structures arise and the condensate becomes unstable as seen 
experimentally |7J. A theoretical description based on the time-dependent Gross- 
Pitaevskii equation (zero-range two-body interaction) was given in jHl E] ■ By definition 
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then only average properties are incorporated although the dominating decay mechanism 
is very sensitive to correlations. 

Correlations are clearly fully included in exact solutions of the full problem as 
attempted in few-body physics. For fermionic systems the current limit is about 10 
particles, see e.g., ^Oj- The simplifications for identical bosons allow computations of a 
larger number of particles especially when using variational methods like the quantum 
Monte Carlo method, see e.g., [TT]. However, detailed investigations for relatively few 
boson or fermion systems already require a substantial effort ^21 El- A larger number of 
particles can be handled if only specific properties are wanted and not the full correlated 
wave function. In particular quantum Monte Carlo calculations have reproduced density 
profiles in agreement with the Gross-Pitaevskii results [T^l I15j. 

The need to account for correlations seems unavoidable in experiments where cluster 
structure is important, e.g., formation of molecular dimer states. The simplest example 
is probably the three-body recombination process, where two of the atoms in the many- 
body system react and form a two-body bound state. In [THj is suggested to apply a 
Feshbach resonance to create a hybrid atomic-molecular Bose-Einstein condensate. The 
atom-molecule coupling is included on top of the usual mean-field equations. It is then 
conjectured, that the ground state of a Feshbach resonant Bose-Einstein condensate in 
reality is a mixed condensate of atoms and di-atomic molecules. Adhikari |T7j studied 
the coupled system of atoms and molecules from the Gross-Pitaevskii equation and 
predicted oscillation phenomena. The experiment reported in [fSj used the tuning of 
a Feshbach resonance with a collapse-burst process as a result. It was subsequently 
shown [T5J E| that the coherent burst-remnant oscillations could be accounted for 
by the presence of a molecular Bose-Einstein condensate. Furthermore, the creation of a 
di-atomic molecular condensate of fermionic atoms was recently observed [221 I2H] giving 
additional evidence for the creation of a mixed condensate. 

Descriptions of correlations within iV-body systems suggest the use of few-body 
techniques which are tools to understand the few-body structures essentially decoupled 
from all other degrees of freedom. This suggests to extend the use of suitable three- 
body formulations. A particularly promising set of calculations were reported in J2] 
for an isolated three-body system with total angular momentum zero. They varied the 
scattering length and described systems with any number of bound two-body states. 
Moreover, they studied a range of excited states and concluded that such higher-lying 
condensate-like states do not collapse under the usual conditions when N\a s \/b t > 0.5 
[21], where b t is the trap length of an external harmonic field. 

Generalization of this work to iV-body systems started in where average 

properties of boson systems were investigated with hyperspherical coordinates. In 
[2H1 12Z| this adiabatic hyperspherical method was extended to explicitly include two- 
body correlations in iV-body boson systems. This structure of the wave function is 
beyond the mean-field. The application for the particle number N = 20 was extended 
up to iV = 10 5 |2H1- Scaling properties were deduced as function of scattering length 
and particle number and analytic expressions were derived for the adiabatic potentials. 
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The results from j2Hl EH EE] indicate that the mean-field properties for dilute 
systems are reproduced in addition to the two-body correlations. Better understandings 
of validity conditions and connections between mean-field models and the adiabatic 
hyperspherical expansion method are desirable. Of particular interest are correlations 
for large scattering lengths, which cannot be studied by mean-field methods. 

Before proceeding it may be useful to express our definition of correlations in 
iV-body systems, i.e., as structures indescribable by mean-field wave functions. For 
example, if two particles form a bound state it is possible to formulate a mean-field 
theory for such dimer states. It is also possible to construct a theoretical formulation 
where mean-field wave functions are used for each species in a coupled system of 
single particles and dimers. Even when dimers can separate and two particles combine 
to dimers this could still appropriately be called a mean-field treatment. Examples 
are the Hartree-Fock-Bogoliubov (HFB) formulation in ^H] for atoms and molecular 
(bound or unbound) dimers and the HFB approximation for nucleons and unbound 
pairs of nucleons frequently applied to the nuclear many-body problem [29]. When 
all divergences are removed by renormalization, the restriction of the Hilbert space to 
mean-field wave functions allow computations for any interaction parameter, e.g., large 
two-body scattering length. Obviously, this does not imply that the true many-body 
correlations can be described, only the structures allowed by the wave function. The 
key point is the allowed Hilbert space in a specific formulation. In [19] both atoms and 
molecular dimers are simultaneously allowed. In [2H1 EH EE] not only molecular dimers 
are allowed, but all kinds of diatomic correlated structures with non-zero higher-order 
correlation functions are included. 

The purpose of this article is to discuss both qualitative and quantitative gross 
properties of iV-body boson systems where two-body correlations explicitly are included. 
In section El we briefly summarize the hyperspherical theory for studying correlations. 
The connections to mean-field methods and other descriptions of correlations in many- 
body systems are not previously formulated and we include a general discussion in 
section l2~2l In sectionElwe present the hyperspherical potentials and discuss the analytic 
parametrization of the effective radial potential extracted in [2E]- We derive scaling 
properties and discuss qualitatively a possible scenario for decay and collapse of the 
condensate after sudden changes of the effective two-body interactions. In section HI 
we discuss details of relations to the mean-field, improvements over the mean-field, and 
ranges of validity of the mean-field and the present hyperspherical method. This section 
contains essentially only new results. Finally, we summarize and conclude in section 

2. Theory 

We study N identical bosons of mass m trapped by an external harmonic field of angular 
frequency uj. We assume interaction via a short-range two-body potential V, which may 
depend on the relative spin state. However, we shall see that the interaction basically 
only enters through the relevant scattering length and the formulation is therefore the 
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same as for spinless bosons. The Hamiltonian is then given by 

N -2 i N 
8=1 i<j 

where r*j is the position of particle i, Pi the conjugate momentum, and r^- = \fj — r£| 
is the interparticle distance. The interaction part is independent of the center of mass. 
Both the kinetic energy and the external harmonic field can be separated into parts 
depending on the center of mass and parts only depending on relative coordinates. For 
this we use the relation 

N N 
i=l i<j 

where R = ^nVi/N is the center of mass coordinate. This immediately leads to the 
convenient definition of the hyperradius p 

N N 

1 _ ' 

~~ N 

i<j i=l 

The relative degrees of freedom are first related by iV — 1 Jacobi vectors ff k 
pJU] . We next choose a new set of coordinates, the hyperspherical coordinates, 
to describe the 3iV — 3 relative degrees of freedom. The hyperradius p sets the 
overall length scale for the system, the angles a k determine the N — 2 relations 
r] k = pcos«jv-i cosaAT_ 2 . . • cosa^+i sina fc between the lengths of the Jacobi vectors, 
and 2(N — 1) angles determine the orientations of the Jacobi vectors. All 3N — 4 
hyperangles are collectively denoted by fl j3H E2] ■ 

2.1. Adiabatic hyperspherical method 

The Hamiltonian then separates into a center of mass part, a radial part, and an angular 
part depending respectively on R, p, and Q |S2] 

H = H c . m . + H P + , (4) 



^£4 = 5> 2 -^ 2 - (3) 



2mp 2 ' 



fcn =AjU + ^X>(>'«). (7) 

i<j 

The angular kinetic energy operator A] y _ 1 is given recursively by 

A 2 P 

Al = Til + + -A— . (8) 

k k cos 2 a k sm 2 a k 



d 2 3fc-6- (3A;-2)cos2a fc d 
oaf. sm 2a k oa k 
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where is the angular momentum operator associated with fj^. 

The relative wave function \I/(p, Q) obeys the Schrodinger equation 

(h - H c . m .)*(p 7 n) = ev( p , n) . (io) 

The adiabatic expansion of the wave function is 

oo 

*(p,n) = ^F v (p)^( P ,n) , F v {p) = p- { " N - A)l2 Up) , (ii) 

where <&„ is an eigenfunction of the angular part of the Hamiltonian with an eigenvalue 
h 2 \ u (p)/(2mp 2 ) 

hn$ v {p,n) = \M$ v {p,Q) . (12) 

Neglecting couplings between the different i/-channels yields the radial eigenvalue 
equation for the eigenfunction f v and the energy E v 

h 2 d 2 



2m dp 



,2 



+ U v (p)-E v )f v (j>) = 0, (13) 



2mU v {p) _ K (3iV-4)(3iV-6) p 2 

h? p 2 + 4p 2 + bt ' K ' 

where bt = y/H/ (muj) is the trap length and the adiabatic potential U v is a function 
of the hyperradius. It consists of three terms, i.e., the external field, the generalized 
centrifugal barrier, and the angular average of the interactions and kinetic energies. The 
neglected non-diagonal terms are for large hyperradii less than 1% of the diagonal terms 
for attractive Gaussian potentials. Thus, the center of mass motion is separated out and 
the hyperspherical adiabatic method is promising simply due to small coupling terms. 
The remaining problem is the determination of the angular potential A from the angular 
eigenvalue equation. 

2.2. The wave function 

We have so far not assumed specific structures or restricted the allowed Hilbert space 
for the many-body wave function. At some point we need to make a suitable ansatz 
for the angular wave function <& u {p, ^)- However, first we shall relate to the historically 
successful approaches to describe a many-body wave function. 

2.2.1. The Hartree mean- field description. The ground-state Hartree product of single- 
particle amplitudes [33J 

N 

*H(rL,f 2 , . . . ,r N ) = JJ^s.p.(^) > (15) 

i=i 

is for the non-interacting gas in the external field given by the amplitude 

iM*) = Ce-^t) , C- 1 = . (16) 
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With the relation J^lLi rf = p 2 + NR 2 this can be rewritten as 

N 

r 2 , • • • , r N ) = C N exp [ - ]T r 2 /(2b 2 ) 

i=l 

= C N e -?/W) e -NiP/W) = To (R)F (p)% . (17) 



The separation of the center of mass motion assures that the ground-state center of mass 
function always is T (.R) = CN 3 / 4 exp[-NR 2 /(2b 2 )}. Then equation (dZJ) is a product 
of the lowest solution for the motion of the center of mass in a trap, and the lowest 
hyperspherical wave function F $o in equation (JUJ), where F (p) oc exp[— p 2 /{2b 2 )] and 
the angular part $o(P; ®) is a constant. This relation in equation (fTTj) between ordinary 
cartesian and hyperspherical coordinates is valid for any length parameter b t . Therefore 
a mean-field product of identical single-particle Gaussian wave functions is equivalent to 
a hyperradial Gaussian and a constant angular wave function, i.e., with no dependence 
on hyper angles Q. 

In reality the interactions produce correlations and the hyperspherical wave function 
deviates from a hyperradial Gaussian multiplied by a constant hyperangular part. 
Therefore the mean-field Hartree product wave function is not exact. However, a 
measure can be obtained by calculating the single-particle density n, given by 

n(n) = J d 3 f 2 d 3 f 3 ---d 3 f N \^(p,n)T (R)\ 2 , (18) 

which can be compared with the mean-field analogue iV's.p.^i)! 2 - The 3(iV — 1)- 
dimensional integral in equation (j!8)l is very complicated with the full numerical 
hyperspherical solution. To get an idea of the possible structures we assume a constant 
angular part $(p, Q). We expand the hyperradial density distribution on Gaussian 
amplitudes with different length parameters af 

\np)\ 2 = ^ r( sjJ ]n sN^ - p2/a ', (19) 



where J2j c j = 1 assures that F(p) is normalized as J °° dpp 3N 4 |F(p)| 2 = 1. This yields 

NT^ 1 2,n2 n (N-l)a 2 + b 2 

^) = E c .^3 e ^' B ' = ~ ( 20 ) 

which is equivalent to (rf) = j d 3 fi n(fi)rf, since 

«> = > 2 > + <* 2 > = §(i - 1) £ ^? + 1& pd 



and 



J j 



3 
2 
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The mean square distance between the particles can then be obtained from the Gross- 
Pitaevskii, or mean-field, approximation for (rf) by the relation 

^-^fa-W)-^(w> -£!«?). <*> 

For a non-interacting gas in a harmonic external field the energy E is related to 
E oc (rf) by the virial theorem. 

These relations are derived and valid only for Gaussian wave functions. However, 
the true mean-field solution is not strictly a Gaussian, although such an approximation 
is rather efficient as pointed out by Pethick et al. [I] . The above results can be used to 
relate an approximate Gaussian mean-field density distribution to a similar hyperradial 
distribution implicitly, assuming constant angular wave function corresponding to 
uncorrelated structure, see Bohn et al. [25 . 

2.2.2. Faddeev-Yakubovskii description. We seek the effect of correlations and have to 
operate beyond the mean-field. Let us first consider the Faddeev-Yakubovskii techniques 
where the proper asymptotic behaviour of the wave functions directly is taken into 
account |3U • This formulation is well suited when the large distance assymptotics 
is crucial as expected for low-density systems. 

Faddeev J33]. initially studied three-particle systems (N = 3) where one of the two- 
body subsystems is bound, and the other subsystems are unbound. He wrote the wave 
function as $ = 0i2 + 0i3 + 023 with the three terms given by suitable permutations of 

= M^y* 1 ^^ 3 , (24) 

where -R23 = (^2^*2 + ^3^*3) / ("12 + ^3) is the center of mass of the bound subsystem and 
K23 is the conjugate wave vector. A generalization of this three-body wave function is 

N 

<fiij = <f>ij{nj) exp (i ^ ittfk + iKijR^ , $ = ^ faj • ( 25 ) 

k^i,j i<j 

When all relative energies are small, Kij ~ and kk — 0, we obtain 0^- ~ (pij(fij). 

Generalization to an iV-particle system was formulated by Yakubovskh who 
arranged the particles into all possible groups of subsystems and thereby formally was 
able to include the correct large-distance asymptotic behaviour for all cluster divisions 
|35j . The decisive physical properties are related to the division into clusters which for 
N = 3 amounts to three possibilities. The three Faddeev components are related to the 
number of divisions and not the number of particles. For N > 3 the number of cluster 
divisions is much larger than N. For iV particles the wave function is therefore written 
as a sum over possible clusters 

$ = ^ ^(cluster) . (26) 

clusters 

This method is often applied with success in nuclear physics |3fi[ l3*7| EHj- In a 
dilute system two close-lying particles are found much more frequently than any other 
cluster configuration. Then the dominating terms in the general cluster expression in 
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equation (|26|) are the two-body clusters where the remaining particles can be considered 
uncorrelated and described by plane waves or in the mean-field approximation. The wave 
function then reduces to the form 

N 

$(p,n) = 5^(p,n). (27) 

i<j 

2.2.3. Jastrow procedure. The Jastrow variational formulation was designed to account 
for correlations [SHI EDI E]- We will briefly comment on the Jastrow ansatz, since it 
provides a physically transparent reason for writing the wave function as a Faddeev 
sum in the dilute limit. A connection between the Jastrow ansatz for the relative wave 
function [H] 

N 

*fon) = n#W ( 28 ) 

i<j 

and the Faddeev formulation is possible. We write the two-body Jastrow component as 
a mean-field term multiplied by a modification expected to be important only at small 
separation, i.e., (omitting normalization) 

^(ftj) = e -^ /(2OT> t 2 )[l + 0(fy] , 0(f) = for r > r , (29) 

where we introduced the length scale r beyond which deviations due to correlations 
vanish. With equation (JHJ) this leads to the relative wave function 

JV 

*(p, Q) = e~" 2 /^ H [l + 0(7%)] (30) 

N N 

= e-fl^[l + Y,m 3 )+ £ + (31) 

i<j i<j^k<l 

Through equation ([17]) the Gaussian mean-field Hartree-ansatz is obtained for a 
non-interacting system in the harmonic external field. For a homogeneous density 
distribution with b t — > oo the mean-field solution is a constant |42) . For a sufficiently 
dilute system it is unlikely that more than two particles simultaneously are close in 
space, i.e., both r^ < r Q and r k[ < r . Therefore the expansion in equation (|3~Tj) can be 
truncated after the first two terms, i.e., 

N N N 

n [i + <Kry ] -i + E = E Liv-D/2 + 0(r ^ } ] • (32) 

i<j i<j i<j ^ 11 

Redefining the two-body amplitude we end up with a Faddeev-like sum as in 
equation (j2"7l) . 



2.3. Two-body s-wave correlations 

The conclusion from the preceding subsection is that a wave function of the form 

N 

*(p,n)=F(p)J2<f>(P,r ij ) (33) 

i<j 
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incorporates both the mean-field properties through F(p) and the correlations in 
addition to the mean-field through the Faddeev-components <fi. We therefore decompose 
the angular wave function $ in equation (jllj) (omitting indices v), into the symmetric 
expression of Faddeev components (ft 



where the last approximation assumes that only relative s- waves contribute, leaving the 
dependence on the distance r^. Higher partial waves could in principle be included but 
the numerical complications would increase rather dramatically. Thus, we have "only" 
assumed relative s-waves between each pair of particles as appropriate for small relative 
energies and large distances. The capability of this decomposition for large scattering 
length has been demonstrated for N = 3 by an application to the intricate Efimov effect, 
which also arises precisely for small energies and large distances [13*1 144j . 

The Faddeev ansatz in equation (j27j) can be formally established as a generalized 
partial wave expansion in terms of the hyperspherical harmonic kinetic energy 
eigenfunctions. The two-body s-wave simplification then appears as a truncation of this 
expansion to include only the lowest hyperharmonics for the description of the remaining 
N — 2 particles. Since this function is a constant we arrive at equation (|34|) . This s-wave 
assumption emphasizes two-body correlations. The method can be extended to include 
higher-order correlations directly in the form of the wave function, e.g., three-body 
correlations, as suggested by |3T] . 

In conclusion, we believe that the Faddeev ansatz with two-body amplitudes 
accounts for the important two-body correlations when the system is sufficiently dilute, 
and at the same time keeps the mean-field information about motion relative to the 
remaining particles. An extension of this technique would be a feasible, but perhaps 
intricate, approach to study three-body correlations in denser systems and in connection 
with the important process of three-body recombination within iV-body systems. 

3. General properties 

The method outlined above leads to the effective radial potential in equation ()14j) and 
the radial equation IjlUjl . This huge simplification is hiding all the complications and 
the detailed information in the angular eigenvalue computations. The key quantity is 
then the function A, which determines the properties of the radial potential. 

The angular eigenvalue equation (fT2|) can by a variational technique be rewritten 
as a second-order integro-differential equation in the variable ckat-i F° r atomic 

condensates the interaction range is very short compared with the spatial extension of 
the iV-body system. Using this short-range property of the interaction in the integro- 
differential equation simplifies even further to contain at most one-dimensional integrals. 
The validity of our approximations only relies on the small range of the potential, 
whereas the scattering length can be as large as desired. 



A' 



A' 




(34) 
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The general structure of the interaction between neutral atoms is an attraction at 
longer distances arising from mutually induced polarization. At shorter distances the 
Pauli blocking dominates and causes the effective interaction to be repulsive. This is 
often modelled by potentials similar to the van der Waals potential. In the present 
formulation it is possible to use any short-range potential also with a finite repulsion 
at the core, e.g., a sum of two Gaussians with repulsion at the origin and attraction at 
larger distances. However, for the large distances crucial for the condensate's properties, 
only the scattering length is important. We therefore first apply a Gaussian potential 
V{r) = Vq exp(— r 2 /b 2 ) and study dependence on the scattering length a s for a fixed 
range b. It is convenient to measure the strength of the interaction in units of the 
Born-approximation a-Q of the scattering length 

a B = ^2 / d 3 r kl V{r kl ) , (35) 

which for the Gaussian potential is a-Q = v /7rm6 3 Vo/(4/z 2 ). Physical results when 
p ^> VNb are independent of the shape of the potential fI7\ . 

The effective two-body interactions can vary enormously for different systems and 
different experiments. Depending on the strength of the interaction the two particles 
may form a bound state of relatively small radius (nm) compared to the typical size 
(/im) of a Bose-Einstein condensate. For the alkali atoms there are usually several of 
such bound two-body states. Scattering of two atoms at sufficiently low relative energy 
depends only on the two-body s-wave scattering length a s . Large distances appropriate 
for dilute systems can then be expected to be determined almost entirely by a s . At 
higher densities also the effective range of the interaction may be significant. 

We then first solve the angular equation by the method of finite differences [27J. The 
basis points are chosen to catch the rapidly varying parts of the wave function and the 
finite short-range potential. This implies that the points vary strongly with hyperradius 
and particle number. With the angular eigenvalue and wave function we then continue 
to solve the much simpler radial equation where only one adiabatic potential is included. 

3.1. Angular potential 

Two-body interactions are responsible for the properties of the many-body system. In 
our formulation, first the properties of the angular eigenvalues are determined and next 
they enter decisively the effective potentials and the radial equations. Qualitatively the 
results depend on the sign of the scattering length and the number of two-body bound 
states. This is understandable as the atoms in a dilute system at low energy effectively 
interact as in a two-body scattering situation. Higher-order processes seldom occur and 
do not contribute to the properties of the dilute system. 

The study in [2H] included variations of the interaction strength and the number 
of particles. We show in figure ^ the lowest angular potential from equation (fT2|) for 
iV = 100 for Gaussian two-body interactions with various scattering lengths. When 
a s = —b (solid line) the potential has no bound two-body states. The lowest angular 
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10 2 10 3 10 4 10 5 

p/b 

Figure 1. The lowest angular eigenvalues A for N = 100 bosons interacting via a 
Gaussian two-body potential V{r) = V exp(— r 2 /b 2 ) with zero or one bound two-body 
states. The scattering lengths a s /b are indicated on the figure. 

eigenvalue is zero at p = 0, decreases then through a minimum as a function of p 
and continues afterwards to approach zero at large hyperradii as a s /p. Increasing the 
attraction (broken lines) decreases all angular eigenvalues for all p-values. The details at 
smaller hyperradii hardly change with large variations of the scattering length. However, 
at larger distances the approach towards zero is converted into a parabolic divergence 
as soon as the scattering length jumps from negative (dotted) to positive (dot-dashed) 
corresponding to the appearance of a bound two-body state. 

Generally, an attractive finite-range interaction can support a certain number A/b 
of two-body bound states for both positive and negative scattering lengths. Then 
the lowest angular eigenvalues, A , \\, ■ ■ ■ , A^/- B _i, describe these bound two-body states 
within the many-body system at large hyperradii, i.e., they diverge to — oo as seen in 
figure [U The next eigenvalue Xj^ B converges to zero at large distance and corresponds 
to the first "two-body-unbound" mode. Through the derived adiabatic potential this 
mode is responsible for the properties of atomic Bose-Einstein condensation, where no 
clusterization is allowed. 

The detailed numerical analysis in |28J resulted in a parametrization for the 
behaviour of these A-functions. Here we restrict ourselves to attractive two-body 
interactions in two different regimes: i) no bound two-body states and a s < 0, and 
ii) one bound two-body state and a s > 0. For hyperradii exceeding a lower limit po, 
which roughly is at the minimum, i.e., p > p = 0.87A^ 1 / 2 (6/|a s |) 1//3 6 the analytical 
expressions are |2H1 

A, W ,H-|A iW p)|(l + ^) 
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X J i-^i-Wtm) whena s <0, 

|A 4 (JV,p)| + |A,(7V,p)| Wnel1 > U ' 

where Xs is the expectation value of fiQ for the zero-range interaction Vg(r) = 
A-n-h 2 a s S(r) jm in a constant angular wave function Q) =constant, i.e., 

2 r(^=^) a s jv>i 3 /3 Ar7/2 a 



w>> - V; Fi#»y 7 - m^j ' (37) 

Aoo(A^) = - 1.59iV 7/3 , (38) 

A ( 2 ) (p) = E ^^-, ^ (2) = — (39) 
at rn\a s \ z 

The number c approaches unity when the scattering length becomes very large. The 
factor (1 + 0.92iV 7 / 6 &/ p) reflects dependence on potential details like the finite range b of 
the Gaussian two-body interaction. At p ~ iV 7/6 |a s | we have Xs ~ Aoo ~ X^ 2 \ For small 
hyperradii p < po we use for all a s the perturbation result obtained as the expectation 
value of the two-body interaction in a constant angular wave function, i.e., 

A. = A«»)(p).«^, fo r p<p„, (40) 

where we use equation (f4T)|) for small p when |A^(p)| is smaller than the expression 
equation (jSHJ). Then the angular eigenvalue A is defined analytically for all p. These 
expressions describe accurately the results of full numerical computations for any two- 
body interaction as soon as p is larger than p . We should emphasize that the small- 
distance region where A a = A^(p) is sensitive to the specific choice of two-body 
interaction. 

The results of the parametrization in equations (|36|) and (|40|) are illustrated in 
figure El for N = 100 for a larger range of scattering lengths than in figure C] in order 
to show the quality of the parametrization. The pronounced deep minimum at p ~ p 




Figure 2. The angular eigenvalue A, equations and (|40|) . for N = 100 as function 
of p for the different scattering lengths given on the figure in units of the range a s /b. 
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is in the region depending on the two-body potential and reflects only the qualitative 
behaviour of the numerically correct lowest angular eigenvalue. After this strongly 
attractive region at small p the eigenvalues approach zero. As the size of the scattering 
length increases the eigenvalue develops a plateau at a constant value Aoo independent 
of a s . Eventually at large p the eigenvalues vanish as \$ for a s < and diverges to — oo 
when a s > 0. 

When a s < the analytic and the correct eigenvalues both exceed the asymptotic 
zero-range result, i.e., A a > \$ for all hyperradii. This means that the true ground state 
energy is higher than the energy obtained with the zero-range interaction. Thus the 
ground state energy from our model is higher than the mean-field energy. The origin of 
this sequence of energies is that the zero-range interaction inevitably leads to diverging 
energies for smaller distances. The present model avoids this non-physical short-range 
collapse. 

When a s > the interaction is effectively repulsive at large hyperradii and we find 
analogously that an analytical expression in this case for the second angular eigenvalue 
obeys A a < A<$ for all hyperradii, due to the divergence of Aj — > +oo as p — > 0. 
Correspondingly we get energies smaller than the zero-range mean-field result in the 
positive Bohn et al. [25J obtained in this case only energies higher than the 

Gross-Pitaevskii energy. 

3.2. Radial potential 

The parametrization in equation ()36|) leads to an analytic expression for the radial 
potential. We can then also study the properties of the radial potential and derive 
physical quantities like the energy and the root-mean-square separation between bosons. 
In particular the attractive two-body potentials generally give rise to a large number of 
negative-energy many-body states. Using the method described in it is possible to 
estimate the number J\f of such bound states, i.e., 



where >( p > denote the negative part of the radial potential U(p). 

3.2.1. Features of the analytic expression. The radial potential obtained from 
equation (J3l)j) is shown in figure |3] as function of the hyperradius for a series of 
different particle numbers and scattering lengths. The strongly-varying short-distance 
dependence is omitted to allow focus on intermediate and large hyperradii. When an 
intermediate barrier is present the condensate is described as the state of lowest energy 
located in the minimum at large hyperradius. This minimium exists for a s < when 
N\a s \/b t < 0.5 as established in [251127] 

The behaviour at very small hyperradii can be constructed from equation (f4T)J) . 
However, now the central value of the two-body interaction enters explicitly and the 
resulting radial potential therefore depends on the short-distance behaviour of this 




(41) 
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interaction. This model-dependence extends to larger distances where the perturbation 
expression in equation (J4~U|) is invalid. This region of p up to \/Nb is not very 
interesting in the present context and we therefore only crudely connected the analytic 
parametrization in equation (|36|) and the expression in equation ()4Uj) to allow extraction 
of the model-independent result. 

Moving alphabetically in figure H3 from (a)-(f) we first in (a)-(d) maintain the 
particle number N = 6000 while only the scattering length a s varies. From (d)-(f) 
we maintain a s /b = —0.35 and vary N. In (a) the two-body interaction is zero, a s = 0, 




p/b p/b 

Figure 3. Radial potentials with b t /b = 1442 and (a) N = 6000, a s ~ 0; (b) N = 6000, 
a s /b = -0.05; (c) N = 6000, a s /b = -0.18; (d) N = 6000, a s /b = -0.35; (e) N = 3000, 
a s /b = -0.35; (f) N = 500, a s /b = -0.35. The dashed lines are with a s = 0. The 
divergence U(p) — > +oo when p — > is not shown. 

leading to a vanishing angular eigenvalue, A = 0. The effective radial potential then 
consists only of centrifugal barrier and external field with one minimum. In (b) we turn 
on an attractive potential, a s = —0.056, sufficiently strong to overcompensate for the 
centrifugal repulsion and create a second minimum in the radial potential at smaller 
hyperradius. An intermediate barrier is left between the two minima at small and large 
hyperradii. A further increase of the attraction in (c) removes the barrier while leaving 
a smaller, flat region. The negative-potential region around the minimum at small 
hyperradius is now even more pronounced. This tendency is continued in (d) with a 
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stronger attraction. 

With the scattering length from (d), a s = —0.356, and a decreasing number of 
particles the intermediate barrier is slowly restored while moving to smaller hyperradii. 
In (e) for iV = 3000 a barrier is about to occur, and in (f) with only iV = 500 an 
intermediate barrier is again present between a minimum at small and large hyperradii. 



3.2.2. Large scattering length A large scattering length implies through equation (J3H|) 
an intermediate region in hyperradius where the angular potential is almost constant. 
More specifically, when p < N 7 ^ 6 \a s \ and A ~ Aoo two of the terms in the radial 
equation ()14p add to a negative value 

A^ + (3iV-4)(3iV-6) <0| 

which implies that no repulsive barrier is present. Then the effective potential is p~ 2 
until the trap begins to dominate. 

We show in figure|3]the analytic radial potential corresponding to one of the angular 
eigenvalues from figured We observe that the radial potential is negative in a large range 
of hyperradii which can be divided into three different regions. For small hyperradii the 
radial potential has a minimum. For intermediate hyperradii the angular potential is 
a constant and therefore the radial potential behaves as — 1/p 2 . This is from figure 121 
seen to appear for p/b between 10 2 and 10 5 . When p > N 7 ^ 6 \a s \ the angular potential 
vanishes as —1/p, so the radial potential vanishes as —1/p 3 , although not clear on the 
figure. Finally the trap oc p 2 dominates with positive contributions at large hyperradii. 
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Figure 4. Analytic radial potential obtained from equations (|14|) and (|36|l for N = 100 
and a s /b = — 10 4 . 



The bound states in this potential can be divided into groups according to their 
hyperradial extension. The total number of such states is written as M = M\ +Me +Ni 
where Mi, Me, and A2 are the number of states located respectively in the attractive 
pocket at small hyperradii, in the intermediate —1/p 2 region and at hyperradii large 
compared with the scattering length. 
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With the analytic expressions inserted into equation ("4~*[j) to obtain the effective 
potential we find a crude estimate of the number of self-bound states in the pocket to 
be Mi ~ 1.3iV 3 / 2 . This number depends on the properties of the two-body potential but 
the iV 3//2 scaling remains unchanged for all short-range interactions. The outer region 
also supports bound states when the trap length b t is sufficiently large, i.e., b t ^> iV|a s |, 
where we analogously find that Mz — 0.78iV 7//6 . This number may be severely influenced 
by the confinement of the external trap, but again the N 7 ^ scaling remains unchanged. 

The intermediate region is only present when the scattering length is relatively 
large. This region exist when 

&«-^g«|a,|, (43) 

where these limits correspond to values of the hyperradius larger than p m ; n = iV 7//6 & and 
smaller than p max = N 7 ' 6 \a s \. The number of states Me located in this region is again 
obtained from equations (JHJ), an d P^jL see (US 120] • This gives 

Ar E ~Mlnf^), (44) 

7T V p min J 

e 5 - A. - ^-^ -") -!« 1.59^, (45) 



where we used equation (j38|) and assume N is large. The number of these bound states 
is then 

Me ^ 0.40W 7 / 6 In (^pj . (46) 

They are located in the region, where the radial potential behaves as 1/p 2 , which is the 
generic form of the potentials giving rise to the Efimov states in three-body systems 
[4*3*1 14*5~| . These states have characteristic scaling properties relating neighboring values 
of energies and mean square radii. The number of states depends logarithmically on the 
size of the l/p 2 -region as in equation (jSJ), analogous to the three-body Efimov states. 
These states were therefore denoted many-body Efimov states [*""""| . 

This estimate assumed that the external trap has no influence on the hyperradial 
potential for p < p max - However, when the trap length b t is sufficiently small, i.e., 
when ptrap = \/Nbt < N 7 ^ 6 \a s \, the extension of the plateau is truncated at large 
hyperradii. The number of states can then be estimated by substituting p max with p trap 
in equation (jUj) . This yields 



(47) 



The plateau can not exist when the external potential dominates already at small 
distances, i.e., for short trap lengths when p trap < N 7 ^b or equivalently iV > iV max = 
(& t /&) 3 / 2 . This maximum number of particles allowing a plateau and the resulting 
Efimov-like states is for a realistic ratio of b t /b = 1442 therefore obtained to be 
iV max ~ 55000. This estimate is rather uncertain but it illustrates that too many 
particles not only exclude stability of the condensate but also the existence of the 
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spatially extended Efimov-like structures which otherwise might play a role in the 
recombination processes. 

The number of Efimov-like states A/e increases strongly with N as seen in table Q 
for b t /b = 1442 for a few particle numbers. These estimates are more precise than in 
previous work where we obtained JVe = 28 for N = 20 for the same parameters |26j . 



N 


10 


20 


100 


1000 


10000 


Ae 


35 


72 


380 


3632 


24810 



Table 1. The number A/e of Efimov-like states for bt/b — 1442 and \a s \ — > oo. 



The energies and mean square radii of the Efimov-like states are related by the 
expressions (IE] 

E " = -2^)J (1+ « 2) - * = < 48 > 

where the exponential dependence on both the strength £ of the effective potential and 
the number of excited states is highlighted. 

Let us assume that the trap length is large and not responsible for terminating the 
plateau at large distance. We can then crudely assume that the mean square hyperradii 
of the first and last Efimov-like states are given by p^ in = N 7 ^ 3 b 2 and p^ ax = A^ 7 / 3 |a s | 2 , 
respectively. Using equation we then obtain the energies of the first and last 
Efimov-like states 

h 2 h 2 

£first ^ ' E ^ ~ ~ o i 12 • ( 49 ) 

2mo z 2m\a s \ z 

which turn out to be independent of the particle number N. These energies remind 
of the kinetic energy scale of strongly bound two-body states and the expression for a 
weakly bound or a resonance two-body energy, respectively. This does not mean that 
the average distance f between two particles in these many-body states also are given by 
b and a s . In fact, f contains an additional iV-dependent factor, i.e., f rs N 2 ^ 3 b, N 2 ^ 3 \a s \ 
for the two cases. These constant energy limits imply that the density of Efimov-like 
states increases with particle number precisely as the interval scales, i.e., iV 7 / 6 . 



3.2.3. Decay and collapse. The Bose-Einstein condensate is intrinsically unstable and 
decays spontaneously, e.g., into lower lying dimer states. Recombination of two particles 
into a lower-lying (bound) state is possible by emission of a photon, but the rate is 
strongly enhanced when a third particle is involved instead of the photon. This three- 
body recombination process has been suggested to be important in Bose condensates 
The outcome of dimers can not be distinguished directly from these very similar 
processes. Molecular formation from two cold atoms, enhanced by tuning the Feshbach 
resonance, corresponds to absorption of a photon and creation of a meta-stable structure. 
The related change of the surrounding medium could lead to instability collectively 
involving many particles, and much faster decays better described as a collapse jH]. 
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In the experiments by Donley et al. [7] a condensate was first created with effectively 
zero interaction, i.e., zero scattering length as in figure Et- The radial wave function is 
then located at relatively large distances in the minimum created by the compromise 
between centrifugal barrier and external field. The attractive pocket at small distances 
is not present and the condensate appears as the ground state in this potential. Both the 
radial potential and the wave function are shown in figure 03 The effective interaction 
was then suddenly changed by tuning a Feshbach resonance [3B] to obtain a large and 
negative scattering length [Zj. The measurement showed a burst and a remnant of 
coherent atoms. This was interpreted and explained as formation of dimers via the two- 
body resonance, a burst of dissociating dimers and a remnant of an oscillating mixture 
of coherent atoms and coherent molecules ^3 1201 EI] • 

In our formulation the effective potential is suddenly altered by a change of the 
underlying two-body interaction. The corresponding new radial potential, shown in 
figure has a pronounced attractive region able to support a number of bound states. 
The initial wave function is no longer a stationary state in the new potential and a 
motion is started towards smaller hyperradii, where it would be reflected from the wall 
of the centrifugal barrier. The system would then oscillate between the centrifugal 
barrier and the wall of the external field. 

This makes the unrealistic assumption that no other degrees of freedom are 
exploited, e.g., the angular dependence of the wave function or molecular bound states 
described by other adiabatic potentials. Thus direct population of two-body bound 
states and resonances are not allowed. This requires in addition inclusion of the adiabatic 
potential asymptotically describing these states. This is entirely possible within our 
model, but it constitutes a major new numerical investigation where coherent atoms 
and molecules, oscillations between them, and (three-body) recombination are studied 
in the same framework. 

In the present work we confine ourselves to the scenario of macroscopic contraction, 
where the density rapidly increases and dimers quickly are produced and subsequently 
ejected from the trap. We shall in the following make qualitative estimates of the three- 
body recombination rate producing the dimers. Whether this process is significant or 
not remains to be seen. 

To study the process we maintain the chosen degrees of freedom described by one 
adiabatic potential. We expand the initial wave function on the eigenfunctions in the 
new adiabatic potential. The dominating states in this expansion are the highest-lying 
Efimov-like states now present because of the large scattering length which produces the 
plateau region and the p~ 2 potential. These states have a similar large spatial extension 
as the initial wave function. The resulting non- stationary wave function provides a 
specific oscillation time. After a quarter of a period the extension of the system has 
reached its minimum. The wave function at this time T is also shown in figure 

The recombination probability increases with decreasing hyperradius due to the 
higher density, i.e., several particles are close in space and therefore much more likely 
recombine into molecular states. The time scale for three-body recombination is given 
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Figure 5. Wave functions / and effective hyperradial potentials U in dimensionless 
units as function of hyperradius for N — 20. The scattering length is zero up to the 
time t = and then suddenly changes to be large and negative at later times t > 0. 
Potentials and the corresponding wave functions are sketched for t — and at a time 
T after a quarter of a period. The horizontal lines show the stationary negative-energy 
states for t > 0. 




by N(t) = N(0) exp(— i/T rec ), where N is the number of atoms in the condensate. This 
function of the average hyperradius p estimated by 



(51) 



Trec = H\a s \*N* • (50) 

This recombination time for the highest-lying Efimov-like states (p ~ N 7 ^ 3 \a s \) can 
then be compared to the time scale for motion in the condensate which is given by 
T trap w 2-k/uj. We find 

T rcc ^ N 2 fN\a s 

^trap 7T \ b t 

Thus, close to the limit of stability established as N\a s \/b t ~ 0.5, we have T rec 3> T trap 
(N ^> 1) and the recombination process is rather slow for these highest-lying Efimov- 
like states. Still the lifetime must in all cases be shorter than for the initially created 
condensate because the density is larger. 

If these Efimov-like states are populated in experiments where the potential 
suddenly is changed from figure to figure EJi,e they could possibly be indirectly 
observed. A signature of this many-body Efimov effect would be observation of the 
diatomic molecules formed in the recombination process and with the estimated rate 
T rec from equation (|51j). The rate should then be inversely proportional to the square of 
the scattering length reached after changing the potential. The dimers themselves can 
probably not be distinguished from this and other processes, but the measured rate of 
dimers can possibly be separated into different characteristic components. 

These Efimov-like states may exist as quasistationary states essentially decoupled 



Structure of boson systems beyond the mean-field 



20 



from all the other many-body degrees of freedom. The recombination time or the 
corresponding width r Efimov = h/T rcc of the Efimov-like states indicates the degree 
of decoupling. Using equations (jIBJ) and (f4T)J) we obtain 

rEfimov 1 / b_\ rEfimov 1 /_~\ 

E n - E~n~x ~ 4vr 2 iV 11 /6 \JJ ' En - E n ^ ~ 4vr 2 iV 11 /6 1 J 

for the first and last Efimov-like states, respectively. The couplings compared to the 
level spacings are small and decreasing with N. Thus the identities of these states 
could be very well preserved within the many-body system. Still their lifetimes due to 
recombination processes can be very large compared to the time scale defined by the 
external field. These negative-energy self-bound many-body states should essentially 
maintain their spatial extension after the external field is switched off. This is in clear 
contrast to positive energy states where only the ion trap prevents expansion. Thus a 
relatively slow time evolution of the density distribution without external field should 
be characteristic for these very weakly coupled many-body Efimov states. 



4. Connections to the mean-field approximation 

The mean-field is often used to describe a condensate. A Hartree product of single- 
particle wave functions describes successfully a Bose-Einstein condensate of a dilute, 
weakly interacting gas of pointlike particles with n|a s | 3 <C 1, where n is the density. The 
mean-field validity condition is then fulfilled, i.e., the mean free path is long compared 
to the interaction range of the system defined by the scattering length. The low-energy 
scattering properties expressed by the scattering length are then clearly decisive. In the 
following we first comment on the choice of interaction and second on the differences 
between the mean-field method and the hyperspherical adiabatic method. Finally we 
discuss the conditions of validity. 



4-1. The two-body interaction 

The choice of the interactions should be consistent with the Hilbert spaces for the 
different methods. In the mean-field treatment a zero-range interaction is often applied 

V s (r) = ^^5(r) , (53) 
m 

where a s is the two-body s-wave scattering length. This limit can be obtained from a 
finite-range potential where the range approaches zero and the strength is appropriately 
adjusted. We use a finite-range Gaussian interaction 

V G (r) = V e- r2 ' b2 , (54) 

where the Born-approximation as to the scattering length then is a measure of the 
strength, see equation ([35)1 . The Gaussian is in the limit when b — > a representation 
Ab(f) of the Dirac 5-function 

A b (f) = -— y— — e -r2//fe2 , 1=1 d 3 fA b (r) . (55) 



7T 



3/2^3 
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We rewrite equation ()54|) as 



V G (f) = 7r 3 /W A 6 (rO 




(56) 



m 



which has the same form as equation (|53|). but with a-Q instead of a s . Then for a s = a-Q 
we have 



However, a s = ob is only valid when \a-B\/b 0, which is rarely the case. 

The limit of vanishing range b can be reached in several ways, e.g., as in 
equation ()56)1 with a constant a-Q or with an adjustment of Vq to keep a constant a s . 
These limits differ enormously and the optimum choice depends on the purpose and the 
Hilbert space restricting the wave function. If the low-energy scattering properties 
are crucial the constant a s seems to be the choice. However, this does not lead 
to equation (jB^j) . but the strength of the interaction should instead approach zero 
linearly with b. In fact, the scattering length is not even defined for the interaction 
in equation (J53|) . Still, the aim of computing reliable energies in the mean-field 
approximation can be achieved with this strength for dilute systems jlHj. The interaction 
and the Hilbert space must be consistent, i.e., a renormalized interaction follows a 
restricted space to produce the correct energy. In this case the Hilbert space is restricted 
to the mean-field product wave functions. Any extension to include features outside this 
restricted space, for example two-body cluster structures, would be disastrous [SU]. In 
other words the wave functions are very difficult to improve even in perturbation theory. 
Maintaining the finite-range interaction with the correct scattering length then results 
in different properties of the interaction even when the range approaches zero on any 
scale defined by the physics of the problem. Thus the mean-field product wave function 
with a realistic two-body potential would also lead to disastrous results. 

Clearly, the full Hilbert space with the correct interaction must produce correct 
results. Whether the realistic interaction combined with our choice of the space including 
two-body correlation amplitudes can reproduce the main features is not apriori obvious. 
However, the investigations summarized in the previous section demonstrate that the 
energy of the mean-field approximation for dilute systems is reproduced and the correct 
large-distance behaviour is at least approximately obtained. This asymptotic behaviour 
is determined by the scattering length which only implicitly is contained in a given 
combination of range and strength of the Gaussian interaction. This implies that our 
Hilbert space must account properly for the crucial correlations necessary for an accurate 
description at large distances. 

4-2. Hyperspherical formulation with the zero-range interaction 

A reformulation of the mean-field in hyperspherical coordinates was given by Bohn 
et al. [25J. They assumed an angular wave function, where all correlations are 
neglected, and a S- interaction, equation (|H3j) . is used precisely as in the mean-field 



lim Vq{t) = Vs(f) . 



(57) 
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approximation. This results in an angular potential produced by the angular eigenvalue 
Xs in equation (p?7j). With this hyperspherical potential they solve the radial equations. 

Roughly speaking, our angular potential arising from the Gaussian interaction is 
above Xs when a s < 0, and below when a s > 0. The "exaggeration" in of the 
zero-range interaction is a result of including the a s /p-divergence of Xs also for small 
distances. When p approaches zero or the scattering length diverges, these and other 
mean-field methods yield disastrous results. 

The mean-field interaction energy can be estimated as the expectation value of 
the 5-function interaction in equation with a Hartree wave function of Gaussian 
single-particle factors: 



1 



m = ^w e ' r,mv • (58) 

Em = mzll f d% WMrMrU = 2N{N ~ if" , (59) 
2 J y/irmbt 

where we used b t as the size parameter for the wave function since the confinement 
is due to the trap. This wave function is then the lowest harmonic oscillator solution 
obtained without any two-body interaction. 

With hyperspherical coordinates this interaction energy is then related to the 
angular eigenvalue: 

h 2 x 5 (p ) 

2mp 2 

where / is the normalized radial Gaussian function corresponding to the Hartree form 



Bin, = / d P r(p)^g-M , (60) 



M = Ji^pp f >3N - 1V2 ■ (61) 

This radial wave function is not the correct solution obtained by using the effective 
potential corresponding to Xs- However, only this Gaussian approximation allows 
an analytic comparison between the hyperspherical and cartesian mean-field wave 
functions. 



4-3. Properties of the wave functions 

The Hartree wave function is closely related to the hyperradial function in the dilute 
limit and the Jastrow correlated wave function is closely related to the Faddeev-like 
decomposition of the wave function. A direct comparison of the wave functions is not 
possible in general as this requires an expansion on a complete set of basis functions 
in one of the coordinate systems. The necessary calculations involve non-reducible 3N- 
dimensional integrals. 

Instead we use the indirect relations provided in section I2.2.H where energy and 
average distance between particles are characteristic features of the wave function. For 
a given scattering length the energy E is numerically obtained for a Bose-Einstein 
condensate as a function of the particle number. We then calculate the interaction 
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energy defined as E — E where E = 3Nhu/2 is the energy of the non-interacting, 
trapped gas. The results are shown in figure |H1 For attractive potentials the mean-field 
has a local minimum at large average distance and much lower (diverging for zero-range) 
energies at small average distances. The mean-field (quasistable) solution is located in 
the minimum at large average distance. This minimum becomes unstable for sufficiently 
large particle numbers. In the example of figure no stable mean-field solution exists 
for iV > 1000. This is consistent with the stability criterion of about iV|a s |/6 t < 0.55 as 
seen from the x-axis exhibited at the top of the figure. 
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Figure 6. a) Gross-Pitaevskii energy as a function of N for a s /b ~ —0.84 and 
bt/b = 1442. Also shown are the hyperspherical calculation for three particle numbers 
and a s /b = —0.84 (a^/b — —0.5). The dashed line shows the Gross-Pitaevskii energy 
for a s /b = —0.5. The ./V|a s |/&t-axis above only applies to a s /b = —0.84. 
b) Mean square distance between the particles for the cases of a). 



In the same figure we compare to results obtained with the present method for three 
different particle numbers. The interaction energies are remarkably similar to those of 
the stable mean-field solution where the scattering length in the Born approximation 
equals the correct value. We also show the results of the less attractive zero-range 
interaction where the scattering length in the Born approximation is the same as for 
the finite-range potential. Now the mean-field interaction energy is much less negative. 
We should emphasize that this comparison does not include the negative-energy states 
supported by the attractive pocket at short distance, see figure EJ They would appear 
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below the "condensate-like" state shown in figure mo- 
using equations (jUJ) and (J23j) we compare in figure Et> ( r ?2) f° r the solutions 
of the mean-field approximation and the hyperspherical methods. The mean square 
distance decreases with increasing particle number for all calculations with an attractive 
potential. As N approaches 1000 the Gross-Pitaevskii mean-field radius approaches zero 
due to the unavoidable collapse. The same behaviour is seen for radii and interaction 
energies, i.e., the average distance between particles decreases until the condensate 
collapses and the size vanishes in the mean-field while many-body bound states with 
smaller extension play a role in the hyperspherical description. Then also higher-order 
correlations can be expected to be essential and result in recombination processes. 

The average distance is related to the interaction energy E — E . In a harmonic 
trap the relation Eq oc (r 2 ) is valid. For condensates the trap determines the average 
properties. It is then not very surprising that the numerical calculations of (r 2 ) 
show that the interaction energy roughly is proportional to the mean square radial 
difference between interacting and non-interacting systems, i.e., E — Eq oc (r 2 ) — (r 2 )o- 
The similarity of these two sets of second moments indicates that the corresponding 
wave functions also are similar. For weak interactions (very small scattering lengths) 
a stationary many-body state can be approximated by a product of single-particle 
amplitudes. However, stronger attraction between particles must invoke other degrees 
of freedom like clusterization. Then a simple single-particle description is not valid. 

4-4- Validity conditions for the models 

Validity criteria for our model and the mean-field approximation, both for zero and 
finite-range interactions, can be compared for a Bose-Einstein condensate where the 
wave function is located at hyperradii p ~ VNb t . Accurate angular eigenvalues in this 
region are crucial for a proper description. If these hyperradii are sufficiently large, i.e., 
p ~ \fNb t > N 7 ^\a s \, the angular eigenvalue has reached its asymptotic value where 
A ~ A^. This condition is equivalent to N\a s \/b t < iV 1 / 3 which is obeyed by stable 
condensates where iV|a fl |/6 t < 0.5 < iV 1/3 

The different models are valid if appropriately designed, i.e., our model should 
reproduce the correct scattering length, whereas both the zero and finite-range mean- 
field interactions should reproduce this same correct scattering length but by using the 
Born approximation. The interaction energies and sizes would all be similar for the states 
corresponding to the condensate. To make this comparison and reach this conclusion we 
have to assume that the angular wave function is a constant and that the hyperradial 
function is equivalent to the single-particle product in mean-field computations, see 
section l2~2l Otherwise the direct connection between wave functions and their properties 
is impossible. This assumption about a specific form of the angular wave function is 
similar to that of spherical Hartree-Fock computations for identical fermions. 

If we for a given average hyperradius p, through equation (jH}, relate the mean-field 
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average distance r by r pa pj yN, then the density n of the system is given by 

3 3iV 3 / 2 , . 

n pa — — ^ — — . 62 
Anr 3 Anp 3 

The zero-range mean-field method is usually claimed to be valid for condensates when 
47rn|a s | 3 /3 <C 1, see |421 15] . Then the number of particles within a scattering volume 
47r|a s | 3 /3 is on average much smaller than one. 

On the other hand, in the zero-range asymptotic region of p > N 7 / 6 \a s \ we have 
np 3 > iV 7 / 2 72|a s | 3 immediately implying that n|a s | 3 < 1/N 2 <C 1, which means that 
the system is very dilute and both zero and finite-range mean-field energy is accurate. 
For p < N 7 / & \a s \ the large-distance asymptotics are not valid and the zero-range 
mean-field description breaks down. For N l l 2 \a s \ <C p < N 7 ^\a s \ or equivalently 
1/N 2 < n|a s | 3 <C 1 the finite-range, but not the zero-range, mean-field is valid. For 
even smaller distances of p < iV 1//2 |a s | also finite-range mean-field becomes invalid. 

The present adiabatic hyperspherical method with two-body correlations explicitly 
allowed in the form of the wave function is first of all valid in the same region as 
the finite-range mean- field approximation, i.e., for N l l 2 \a s \ < p, where correlations are 
expected to be insignificant. However, the validity range of the hyperspherical method 
with two-body correlations incorporated extends to hyperradii smaller than iV 1 ' 2 )^], 
where two-body correlations are sufficient to describe the clusterizations. 

When higher-order clusterizations occur, any method without correlations higher 
than two-body breaks down. The density when this happens for this hyperspherical 
method is not easily derived. The lower limit is probably when the distance between 
two particles on average equals the interaction range b, i.e., N 1 ^ 2 b < p. However, for 
nuclei with identical fermions the radius at saturation is about iV 1//3 6 where the mean- 
field approximation is very successful. This limit would then correspond to N 5 ^ 6 b < p, 
but identical boson systems may allow even smaller hyperradii. 

In conclusion, the validity regions for the two-body correlated method 
(hyperspherical), finite-range methods (finite-range mean-field and hyperspherical), and 
the zero-range mean-field are estimated to be 

p > VNb for two-body correlated method, (63) 
p > y/~N\a s \ for finite-range methods, (64) 
p > N 7 ^ 6 \a s \ for finite- and zero-range methods. (65) 

These relations can with equation (}62|) be expressed via the density 

n|a s | 3 < f — 7— ) for two-body correlated method, (66) 



n|a s | 3 < 1 for finite-range methods, (67) 

n|a s | 3 < — for finite- and zero-range methods. (68) 

When the density is low, the three approximations are valid and the energies are 
similar. This assumes that the renormalization is appropriate. For higher densities the 
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importance of correlations increases and the mean-field approximations break down. At 
even higher density also two-body correlations are inadequate and the particles want 
to exploit higher-order correlations. In any case, the wave functions can not be better 
than the Hilbert space they span, no matter how precise the energy is computed. 

5. Summary and conclusion 

The method of hyperspherical adiabatic expansion is briefly sketched for a system of 
identical bosons. The form of the wave function is chosen as the s-waves in a partial wave 
expansion of the Faddeev-Yakubovskii cluster amplitudes. This restriction is expected 
to be accurate for large distances and dilute systems. We relate to the Jastrow ansatz 
designed to deal with correlations in rather dense systems. We discuss the theoretical 
connections between these approaches and the mean-field approximation both with zero 
and finite-range interactions. 

The angular eigenvalues in the hyperspherical adiabatic expansion appear as 
crucial ingredients in the radial potentials. We use the analytic expressions recently 
parametrized to reproduce the results of full numerical computations. We first 
discuss the general properties of these eigenvalues as functions of hyperradius for 
arbitrary particle number and arbitrary scattering length. The large- distance behaviour 
corresponding to the zero-range mean-field result is obtained. 

The radial potential has a minimum at large distance when particle number times 
scattering length divided by trap length is less than about 0.5. The wave function of 
the condensate is located in this minimum. In addition, for sufficiently large scattering 
lengths an intermediate region appears with a radial potential decreasing inversely 
proportional to the square of the hyperradius. This region supports the many-body 
Efimov-like states. At much smaller distances a pronounced attractive pocket is present 
when the two-body potential is attractive. We give analytical estimates of the number 
of bound states located in these different regions. We then discuss the decay properties 
eventually arising from recombination processes. In particular the highest-lying Efimov- 
like states located at large distances recombine corresponding to widths much smaller 
than the level spacing. These peculiar states could then leave observable traces. 

Finally, we discussed the connection between this work and the mean-field approx- 
imation. We first emphasized that the effective two-body interactions must be related 
to the Hilbert space for the wave function. We specify the necessary renormalization for 
the mean-field restriction. Numerical comparison for energies and radii are then pre- 
sented. The validity conditions for the models are discussed and expressed as regions in 
hyperradius. These regions increase from zero via finite-range mean-field approximation 
to the hyperspherical adiabatic expansion method. Most of the results are independent 
of the structure of the two-body interaction. The conclusions are derived in terms of 
scattering length, number of particles, external field frequency and occasionally the ef- 
fective range of the two-body potential. 
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